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A  mixed  layer  dynamic  irodel  fcr  t.he  structure  and  evolution 
cf  aeicscls  in  a  marine  regime  is  presented.  The  a=rosol 
spectrum  is  divided  into  continental  and  marine  components, 
with  a  constant  vertical  continental  profile  assumed;  both 
in  and  atove  the  mixed  layer.  The  aerosol  spectrum  is  trans- 
formed into  a  reference  relative  humidity  of  80^.  The 
tempcral  evcluticn  cf  the  aerosol  spectrum  is  predicted  from 
rate  equations  which  require  a  specification  of  the  surface 
production  rate,  the  entrainment  rate  (We  )  »  and  the  irixed 
layer  depth  (h) .  The  model  was  tested  against  the  data  set 
obtained  in  the  Monterey  Bay  during  the  MAGAT  89  experiment. 
The  model  was  initialized  with  both  observed  MAGAT  data  and 
an  equilibrium  initial  value,  gathered  from  the  JASIN 
Experiment.  The  model  was  run  for  radii  equal  to  0.8,  2.0, 
5-0  and  10-0  microns.  The  significance  between  the  observed 
initial  values  and  the  JASIN  data  is  that  the  MAGAT  data 
were  observed  from  the  same  air  mass  as  the  initial  atmos- 
pheric data  and  verification  aerosol  data.  The  JASIN  initial 
aerosol  data  are  based  on  an  equilibrium  state  as  a  function 
of  only  wind  speed  and  reference  relative  humidity.  The 
model  ccntinuously  generated  a  correct  gain  or  loss  of 
aercscl  ccncentra ticns  as  defined  by  the  observed  MAGAT 
data,  and  in  most  cases  the  model  output  is  within  one  crder 
cf   maanitude  of    the    cbserved   values. 
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There  is  an  increased  emphasis  in  ths  ability  to  fore- 
cast the  behavior  and  evolution  of  aerosols.  The  scattering 
and  absorption  of  light  are  influenced  by  aerosol  distribu- 
tions in  the  atmospheric  boundary  layer  and  affect  the 
performance   of      optically   guided      weapon   systems.  The    Air 

Force  is  particularly  interested  in  how  aerosol  extinction 
affects  the  use  cf  precision  guided  munitions  (PGM) 
£ottrell  et  al,197jj/  .  The  Department  of  Defense  (EOD)  has 
PGM  *  s  that  cperate  at  differing  wavelengths  which  range  from 
the  visible  to  the  microwave  regions.  PGM's  have  a  greater 
ability  ic  hit  a  target  than  conventional  munitions; 
however,  the  controlling  factor  is  the  ability  of  the 
guidance  system  to  "see"  the  target.  This  ability  is  depen- 
dent en  the  wavelength  for  which  the  sensors  are  designed 
and  the  properties  in  the  intervening  atmosphere.  The 
degrading  properties  in  the  atmosphere  are  principally 
Eolecular  absorption  and  aerosol  scattering.  The  wavelengths 
for  the  guidance  systems  are  designed  such  that  the  mclec- 
ular  absorption  is  minimized:  therefore,  scattering  ty  aero- 
sols becomes  the  main  concern  once  a  suitable 
absor pticn-free    window   has    been   selected. 

The  model  under  consideration  includes  the  behavior  of 
larine  aerosol  components  as  well  as  previous  continental 
components.  Estimating  the  influences  of  marine  aerosols  on 
elect ro-cptical  (EO)  systems  has  been  studied  /Barnhart  and 
Streete,  1970/  .  Particle  sizes  of  interest  are  those  asso- 
ciated with  locally  generated  sea  salts  because  cf  their 
effects  on  IR  as  well  as  visible  wavelengths.  The  size 
distributions  of  sea  salt  particles  show  a  variance  of 
several    crders    of   magnitude. 


The  ability  to  forecast  the  behavior  of  aerosol  extinc- 
tion from  synoptic  scale  patterns  would  help  in  the  deciding 
which  type  cf  weapon  system  to  employ.  Because  some  systems 
are  launched  frcm  the  air,  it  is  important  that  such  a 
profile  include  the  vertical  distributions  of  aerosol  ccrxi- 
cles.  Models  exist  for  estimating  vertical  extinction 
profiles,  but  they  have  not  been  sufficiently  verified.  To 
do  this,  profiles  cf  actual  aerosol  dara  must  be  gathered 
and  ccmpared  with  the  model  forecasts. 

Mcdels  in  current  use  are  based  on  parameter izaticns  of 
the  effects  of  relative  humidity  and  wind  speed  on  the  equi- 
librium aercsol  distributions  /Wells,  at  al,  1977^  .  Hecent 
evaluations  have  shown  that  these  models  are  limited  to  mean 
distr ibuticr.s  (i.e.,  the  average  aerosol  concentration  at  a 
given  wind  speed  and  humidity)  [Fairall  et  al,  1982aJ  . 
Models  are  limited  because  seme  processes  in  the  atmospheric 
nixed  layer  which  affect  aerosol  concentrations  are  not 
considered,  namely  entrainment  and  inversion  height  changes. 

The  purpose  of  this  study  is  to  present  and  evaluate  a 
model  which  includes  The  meteorological  processes  that 
adequately  describe  the  whele  marine  atmospheric  boundary 
layer  (MAB1) .  An  inversion  represents  a  cap  to  the  vertical 
transport  of  surface  generated  aerosols,  and  is  not 
accounted  fcr  in  previous  models.  The  top  of  the  boundary 
layer  is  capped  by  the  marine  inversion,  where  entrainment 
cf  overlying  air  takes  place.  Because  entrainment  mixes 
clear  (ncn-marine)  air  into  the  marine  layer,  this  process 
is  as  important  as  surface  layer  fluxes  in  determining  equi- 
libriun  concentration.  This  entrainment  process  is  included 
in  the  mcdel. 

Evaluation  of  the  ircdel  cutput  will  be  done  with  the  data 
from  an  experiment  entitled   Marine  Aer_c;s_ol   Generation  and 
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Transport  (MAGAT).  The  experiment  was  conducted  in  the 
vicinity  of  Monterey  Eay,  Ca,  during  the  period  of  28  April 
to  9  Eay  1980.  The  purpose  of  the  experiment  was  to  examine 
the  ccmpatibilit y  of  optical  and  micrometeorolcgical  propa- 
gation theory,  and  tc  extend  dynamic  models  of  the  evolving 
MABL        tc        include        aerosol  and        turbulence        profiles 

{Fair all,  1980  and   Fairall  et    al,1980j     •.       Two    platforms,      the 
R/V    ACANIA,    and    an    aircraft    were    used. 

In  this  study,  twc  24-hour  periods  were  chosen  for  =valu- 
ation;  3  May,  beginning  at  1200  PDT  and  6  May,  beginning  at 
1800  EET.  These  periods  were  chosen  because  the  boundary 
layer  was  undisturbed  (no  fronts  closer  than  100  nm)  for  24 
hours  prior  to  the  starting  ximes.  Additionally,  the  times 
were  selected  due  tc  the  proximity  of  the  aircraft  and 
surface  ship  during  the  experiment  times.  The  approach  was 
to  describe  the  synoptic  conditions  from  24  hours  in  advance 
cf  the  mcdel  forecast  through  the  end  of  each  forecast,  and 
compare   the    model   output  with    the   actual    MAGAT   findings. 
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II.     BACKGROUND 

*.       DESCBIF1ION    OF    EIISTIHG    HCDELS 

Current  models  fcr  estimating  aerosol  equilibrium  distri- 
butions use  meteorological  inputs  of  (10-m  height)  wind 
speed  ard  relative  humidity.  These  two  quantities  are 
considered  because  of  their  role  in  generation  and  transport 
(wind)  and  growth  (humidity)  of  aerosols.  The  Navy's 
Sell s-M una- Katz  (WMK)  /wells  st  alf  197^  is  an  example  of 
this    ccncept.  The    performance   of      this   type  of      model   has 

teen  studied  with  data  obtained  in  the  northern  Atlantic  and 
eastern  Pacific  Ocean  areas.  The  model  output  compared  with 
these  data  is  shewn  in  figure  2.1,  depicting  a  height  depen- 
dence cf  tctal  aerosol  volume  from  a  sample  set  of  eastern 
Pacific  data.  These  profiles  correspond  to  (1)  the  cfcserved 
sea  salt  vclume,  V;  (2)  the  observed  sea  salt  volume 
adjusted  to  80%  relative  humidity,  V0  ;  and  (3),  the  WMK 
predicted  volume  adjusted  to  80ft  relative  humidity 
(circles) .  it  is  clear  f rem  this  figure  that  within  the 
mixed  layer  the  observed  decrease  of  aerosol  volume  with 
height  is  less  than  the  model  predicts.  The  surface  gener- 
ated aercscls  appear  to  be  well  mixed  below  the  inversion 
when    normalized   to    80?   relative   humidity. 

The  assertion  that  existing  models  can  predict  only  a 
mean  value  appears  in  the  results  obtained  in  the  North 
Atlantic,  from  the  JASIN  Experiment.  Figure  2.2  compares  a 
single  radius  size  cf  five  microns  and  the  corresponding 
model  prediction.  The  values  and  trends  in  the  predicted  and 
mean    results      are   in      reasonable   agreement.  However,      the 

standard   deviation    is   three    times   the    mean;    if  one    assumes   a 
normal    distribution,    only   67%   cf   the   observed  aerosol 
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Figure   2.1 


Height    dependence   of  aerosol   volume    and  relative 
humidity  above   the  ocean. 
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Figure  2.2    Aerosol  volume  spectrum  at  r=5  microns 
as  a  function  of  wind  speed,  0. 
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distributions  will  te  withir  a  factor  of  three  of  the 
average.  This  comparison  emphasizes  the  point  that  nc  matter 
how  accurately  a  model  predicts  an  average  aerosol  density 
at  a  given  wind  speed  and  relative  humidity,  the  factor  of 
three  EMS  variations  cannot  te  eliminated  without  consid- 
ering   mere    ireteorolcgical  parameters. 

E.       TEE    SES    BOUNDARY    IAYEB    HODEL 

The  distinguishing  feature  of  this  model  is  the  charac- 
terization of  the  MAEL,  which  is  convectively  mixed  up  to  a 
height  h,  and  capped  by  an  inversion.  The  atmospheric 
profile  is  depicted  in  figure  2.3a,  representing  a  cloud- 
free  mixed  layer  where  water  vapor  mixing  ratio  (g) ,  and 
virtual  potential  temperature  (ev)  are  "well-mixed",  ie, 
independent  of  height  in  the  mixed  layer.  The  assumed 
vertical  aerosol  profile  is  shown  in  figure  2.3b.  The  model 
produces  a  24-hcur  time  evclution  of  an  aerosol  spectrum, 
requiring    a   prediction   of  the    following   at   each   time   step: 

1)  surface   production   rate  of   marine   aerosols 

2)  entrainment  rate   at   the   top   of   the    mixed   layer 

3)  nixed    layer  depth. 

retails  cf  how  these  parameters  are  input  into  the  model  are 
discussed  in  chapter  IV.  The  mathematical  relationship 
between  the  time  rate  of  change  of  the  aerosol  volume  spec- 
trum, dv/dr,  and  these  three  parameters  is  shown  in  the 
following   eguaticn 

dv/dt  =(<wV>0  -  (we  +  wkm)  v)/  h 

1  2  3  (2,1) 

The  model   predicts  the  evolution   of  aerosol  at   five  radii 
(0.8,  2.0,  5.0,  10.0,   and  15.0  microns)   of  both  the  conti- 
nental and  the  marine  (sea  salt)  components  in  the  mixed 
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layer.    The    actual  tine   evoluti.cn    cf   each    radius    is   output   in 
the    fcrm   cf    dv/dr,    with    the    units   of    um  2   /cm3. 

The  model  assumes  r.c  appreciable  concentration  of  marine 
aerosols  atcve  the  mixed  layer.  The  sntrainment  cf  air  into 
the  mixed  layer  will  not  cause  an  elament  of  these  marine 
particles  tc  "escape"  into  the  free  troposphere.  The  ccean 
acts  as  a  source  for  marine  aerosols,  primarily  through  the 
generation  cf  white  caps.  This  is  the  only  input  for  the 
marine  component.  The  ocean  also  acts  as  a  sink  for  the 
continental  aerosols,  which  are  generally  of  a  smaller 
radius.  The  entrair.ment  process  at  the  top  of  the  mixed 
layer  irixes  clear  (non-marine)  air  into  the  MABL.  This 
process  could  be  as  important  as  the  surface  layer  fluxes  in 
deteririning   an    aerosol  concentration. 

C.       MODEL    INPUTS    AND    AEROSOL    INITIALIZATION 

The  model  is  designed  so  that  with  the  exception  of 
initializing  the  aerosols,  all  of  the  calculated  parameters 
are  rasec  en  inputs  from  surface  based  observations,  and 
representative    soundings.    The   input   consists   of: 

1)  Ncn-meteorological  inputs  of  latitude,  Julian  day,  and 
the  local  start  time.  These  are  used  for  a  diurnal  radia- 
tive heating/ccol ing  package.  Diabatic  warming  has  an 
obvious  impact  on  the  life  and  strength  of  the  mixed 
layer.  Additionally,  sea-surface  temperature  (SST)  is 
input  at  the  start  time.  Op  to  ten  forecast  SST  values 
can    te   input    during  the   24-hcur   period   as   well. 

2)  Surface  wind  speed  and  direction  at  the  start  time  and, 
as    with   SST,    up  to    ten  additional   forecast    values. 

3)  Frcir  figure  2.3,  temperature  and  relative  humidity  are 
parameterized    by      virtual    potential   temperature      and   water 
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vapcr  nixing  ratio,  respectively.  Mixing  ratios  are  input 
above  the  mixed  layer,  as  well  as  the  lapse  rate  above  the 
"  j  u  up "  . 

An  initialization  cf  the  six  radii  concentrations  in 
three  regimes  is  required.  To  do  this  it  is  necessary  to 
understand  the  actual  distributions  of  the  continental  and 
marine  components  in  the  mixed  layer.  For  the  purposes  of 
evaluating  the  model,  assumptions  must  be  made  en  the 
distributions  of  the  acquired  MAGAT  data.  The  aircraft 
instruments  could  net  distinguish  the  chemical  make-up  of 
the  individual  particulates;  therefore,  it  is  not  kr.cwn,  for 
a  given  iadii  concentration  in  the  mixed  layer,  hew  much  is 
continental  and  how  much  is  marine.  The  concentrations  of 
each  size  were  calculated  in  the  mixed  layer,  and  then  abeve 
the  mixed  layer.  Based  on  the  assumption  that  there  are  no 
uarine  concentrations  above  the  mixed  layer,  we  conclude 
that  the  concentrations  calculated  from  above  are  all  conti- 
nental, and  therefore  linerally  subtracted  from  the  concen- 
trations in  the  nixed  layer,  leaving  only  marine 
concentrations. 

A  further  adjustment  of  the  input  aerosol  values  is 
required,  based  on  changes  in  relative  humidity  durinq  the 
data  collection  time.  We  have  stated  that  the  growth  of 
aerosols  is  a  function  of  relative  humidity.  Therefore,  a 
distribution  of  aeroscls  gathered  at  90%  humidity  cannct  be 
directly  cenpared  to  another  distribution  gathered  at  80  % 
humidity;  a  reference  humidity  is  required.  Consider  the 
aerosol    volume   spectrum 

v(r)  =  4/37rr3n(r)  (2-2> 


1S 


where    n  (r)    is   the   nuirter   density   spectrum.      V  (r)       is    defined 

as   the      volume    of   aerosol   particles      par   cm3    at  a    reference 

saturation    ratio      of    S=80%.       A      humidity    growth  factor    G (S) 
/Fair  all,    et    al,     1982aJ      is    defined   as 

G(s)-.81  exp(.066S/(1.58-S))  (2,3) 


There  are  further  considerations  concerning  advecticn. 
Clearly  all  aerosols  do  not  originate  locally;  both  above 
and  telow  the  inversion  aerosols  are  advected  into  a  local 
region.  In  terms  of  aerosol  density,  entrainment  acts  as  an 
aerosol  flux  out  of  the  boundary  layer  because  the  concen- 
trations above  and  below  the  layer  are  different.  In  the 
model,  the  entrainment  acts  en  the  "jump"  across  the  inver- 
sion. For  the  purpose  of  evaluating  the  the  evolution  of  a 
local  concentration  of  aerosols,  the  model  will  neglect 
hori2cntal  advection  and  further  assume  a  negligible  local 
production  of  the  continental  component.  The  model  takes 
into  consideration  a  Stokes  gravitational  fallout  term,  wk 
feu,  1979/  .  The  fallout  rates  above  and  below  the  inversion 
are  different  because  of  the  change  in  aerosol  spectra, 
caused  by  the  humidity  growth  factor.  The  Stokes  velocity 
/Fairall,  et  al,  1S£2a/  is  calculated  from  the  following 
equation  where  pw  is  the  density  of  the  droplet  and  E  is  the 
kinematic  viscosity  cf  the  air. 

W„  =  2g(pw-p)rJG2(S)/(9Ep)  <2'U) 
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E.   ECtJIUBFIOM  AEROSCL  HODEL 

As  a  means  of  comparison  of  aerosol  initialization 
scheires,  the  model  offers  the  option  of  initialization  with 
equilibrium  values.  These  initial  values  are  based  on  data 
collected  from  the  JASIN  Experiment  /fairall,  et  al,  1982J  . 
A  large  aerosol  data  base  was  collected  during  equilibrium 
conditions,  with  respect  to  wind  speed,  during  12  hour 
periods  and  then  normalized  to  a  reference  relative 
humidity.  Aerosol  spectra  were  than  grouped  into  six 
different  wind  speed  ranges.  Fcr  the  equilibrium  initializa- 
tion scheme  in  the  NFS  model,  a  pre-assigned  value  of  the 
five  radii  are  used  as  initial  aerosol  values,  based  on  the 
wind  speed  at  the  model  start  time.  An  example  of  the  equi- 
librium aerosol  spectra  from  the  JASIN  Experiment  is  shewn 
in  figure  2.4.  Note  that  the  reference  relative  humidity  for 
this  experiment  was  87  percent.  The  graph  is  referenced  here 
only  to  shew  an  equilibrium  type  distribution;  the  initial 
equilibrium  and  MAGAT  aerosol  values  will  be  shown  in 
chapter  IV.  Once  initialized  with  either  equilibrium  data  or 
calculated  data,  the  model  will  run  in  the  same  manner. 


J  10'  - 


>  10' 
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RH*87% 


Figure  2. 4    Aercsol  spectra  as  a  function  of  size 
for  selected  wind  speeds. 
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This  equilibrium  scheme  illustrates  the  fundamental 
difference  between  the  NFS  model,  which  is  a  dynamic  model, 
and  a  steady-state  type  model.  A  representative  model  using 
this  equilibrium  approach  is  the  wmk  model,  mentioned  in  the 
beginning  of  this  chapter.  This  type  of  output  describes  the 
narine  aerosol  distributions  as  a  function  of  wir.d  speed 
(surface  generation),  relative  humidity  (growth  factor)  ,  and 
elevaticr  (vertical  variation  of  aerosols  with  height 
assuming  a  steady  state  vertical  transport  process).  The 
matheiratical  representation  of  this  model  /Fair  ail,  et  al, 
1982a]  is 


n(r)  =  (r/a)*1.62  (c^  +  CsV   ) 


F*exp  \^ 


-z/hnF-8.5(r/a) 


■) 


(2.5) 


where 

r=    the    particle    size   in    microns, 

u=   the    wind   speed, 

v=    0.5    for    u  <    4    m/s, 

v=   u    -3.5    for   u    >  4    m/s , 

f-    l-Uv/60)  3, 

y-    0.364-0.  00293v  1-25    , 

Z=    height    above    sea   surface,    m    , 

h    =    scale   height, m    (800    m    for   Z<1000    m)    , 

a=    0.81exp     (0.  06  6s/ (1 .  0  58-s)  )     , 

S=   H/100,    (H=  relative    humidity   in   percent). 

The    ether  constants  are  given,    based  on    the   value   cf    v. 

v,     m/s  C-,  C2  <5 

v    <    7  350  1000  1.  15 

v    >    7  0  6900  0. 29 

The  equation  described  produces  a  number  density  spectrum, 
n(r).  This  relates  to  the  volume  spectrum  V(r),  as  shewn  in 
equation    2.  2. 
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III.  DATA  ACQUISITION  AND  SYNOPTIC  CONDITIONS 

A.   DATA  ACQUISITION 

The  data  were  obtained  frcm  a  ship  and  aircraft  experi- 
ment, Marine  Aerosol  Generation  and  Transport  (MAGAT)  ,  held 
from  28  April  to  9  May  1980,  in  a  region  30  to  50  miles  cff 
the  ccast  cf  Monterey  Bay,  Ca.  Data  gathering  was  dene  with 
instruments  mounted  en  both  an  aircraft,  and  a  surface  ship, 
the  E/V  ACANIA.  Aeroscl  data  were  obtained  frcm  the  Airborne 
Research  Associates  turbo-charged  Bellanca,  using  a  Particle 
Measuring  System  (PMS)  Axial  Symmetric  Scattering  Aerosol 
Probe  (ASSAP)  particle  counter.  All  measured  data  were 
sampled  every  2.5  secends,  with  a  two-scan  average  cf  every 
five  secends.  The  scans  were  ccllected  in  "ladder"  profiles, 
during  which  the  aircraft  made  measurements  at  a  constant 
altitude  fcr  two  minutes,  climbed  to  a  new  altitude,  and 
repeated  a  rew  measurement  run.  The  instrument  utilized  60 
size  channels  from  2.8  to  14.0  micron  radius.  In  most  cases 
the  ladders  extended  from  near  the  sea  surface  (3  m)  up  up 
through  the  well-mixed  boundary  layer,  to  a  few  steps  above 
the  inversion.  The  elevation  was  generally  up  to  5  kilome- 
ters. A  typical  ladder  profile  contains  10  to  14  steps.  The 
step  heights  were  randomly  chosen,  but  an  attempt  was  made 
to  keep  each  step  height  consistent  between  ladders.  Air  and 
dew-pcint  temperatures  were  also  measured  and  used  tc  calcu- 
late relative  humidity  for  the  correction  factor  mentioned 
in  chapter  two.  The  aircraft  also  flew  ascending  spirals  (in 
the  vicinity  of  the  ladders)  during  which  other  metecrolo- 
gical  parameters  were  ccllected.  This  data  yielded  vertical 
soundings  similar  to  those  prcvided  by  radiosondes. 
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There      tiers    t«o      aerosol    instruments      fitted      on    the      R/V 
ACANIA.         Or.e   was  the  FMS      model   CSAS    (classical    scattsring) 
and    the    ether   was   model    ASAS    (active   scattering),    controlled 
by    a    IMS      data    acquisition    systam    (DAS-32)  with    a    computer 

interface.        The   shipboard      system    measured      aerosols    in      90 
different   size    channels   from   0.9    to    14.0   micron   radius. 

It  was  ncted,  during  the  flybys  over  the  R/V  ACANIA,  that 
the  aircraft  measurements  did  not  agree  with  that  of  the 
ship.  The  size  distributions  from  the  aircraft  measurements 
were  consistently  smaller  than  the  ship  measurements,  for 
radii  greater  than  1.C  micron.  In  addition,  the  differences 
increased  with  radius.  Since  the  shipboard  aerosol  system  is 
newer  and  had  a  wider  size  range,  and  better  sensitivity, 
aircraft  aerosol  data  were  corrected  to  agree  with  the  ship 
aerosol    data    /Fair all,  1980    and   Fairall  et   al,  198_oJ     . 

Profiles  of  virtual  potential  temperature  and  mixing 
ratio  were  obtained  from  three  different  sources.  The  reason 
for  usinc  these  parameters  instead  of  temperature  and  dew 
point,  is  that  mixed  layer  inversions  are  more  easily  iden- 
tified with  the  former  variables.  The  sources  were  the 
spiral  flights  from  the  aircraft,  the  radiosonde  launches 
from    NFS    and    from  the  R/V   ACANIA. 


E.       SINCETIC   DATA 

Surface  and  500  mb  synoptic  charts  and  the  GOES  WEST 
satellite  images  were  used  to  evaluate  the  synoptic  condi- 
tions. Charts  are  from  the  NOAA  weekly  series  of  daily 
weather  maps.  In  addition  to  data  collected  from  the  R/V 
ACANIA  and  the  aircraft,  local  weather  data  were  also  avail- 
able   ficm   the  U.S.    Army's  Fritzsche   Field    weather    facility, 
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Figure   3.1 


Surface  and  500  mb  anaylsis  for  the  western  0- S 
at  0500  PDT  2,  3  and  4  May  1980. 


24 


1615  02HV80  35A-2   006^3 '22031  SB6 


Figure  3.2   GOES  West  satellite  imagery,  0915  PDT  2  May,  1980. 
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1745  Q3MV8Q  35ft-4  00362*19161   UC2 


Figure  3.3   Same  as  Figure  3,2  except  1245  PDT ,  3  May. 
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2345  fMMYflO  35A-4   QQ35T  19161  UC2 
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Figure  3.4   Same  as  Figure  3.2  except  1645  PDT ,  4  May. 
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Figure  3.5   Same  as  Figure  3.1  except  5r  6,    and  7  Bay,  1980 
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Figure  3.6   Same  as  figure  3,2  except  1745  PDT ,    5  Hay. 
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2245  06HY80  35A-4   00342  19171  UC2 


Figure  3,7   Same  as  figure  3.2  except  1645  PDT ,  6  May. 
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Figure   3.8        Same   as  Figure   3.2   except    1745   PDTr    7    May. 
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Fort  Crd,  Ca.  Satellite  photographs  and  NOAA  maps  for  the 
two  different  time  periods  (2-4  May  and  6-7  May)  ars  shewn 
in  figures  3.1  through  3.8.  Each  three-  day  sequence  depicts 
the  24  hcur  period  being  considered;  24  hours  before  and 
after.  Atmospheric  soundings  for  the  two  model  runs  are 
shown  in  chapter  IV,  figures  4.1  and  4.2. 


C.   S1N0ETIC  CONDITIONS 

Several  weak  frontal  systems  passed  through  the  area 
during  the  experiment.  Showers  occurred  during  the  first  and 
last  cays,  associated  with  the  fronts.  Low  cloudiness  and 
fog  cccured  during  the  morning  from  29  April  to  5  May,  with 
fog  returning  again  en  the  last  day. 

At  the  beginning  cf  the  period,  the  area  was  dominated  by 
a  slcwly  eastward  migration  of  a  cut-off  low  at  500  mb.  3y 
early  morning  on  2  May,  the  area  was  under  the  influence  of 
a  weak  ridge.  On  3  Kay  the  area  was  under  divergent  flew  at 
the  upper  level.  An  upper  level  low  had  formed  off  the  Baja 
Ca.  coast  en  4  May,  leaving  the  area  under  the  influence  of 
a  col.  Cn  5  May  the  area  was  between  an  upper-level  trough 
and  ridge,  and  by  6  Kay  the  area  was  on  the  back  side  of  the 
trough.  Eecause  cf  the  deepening  of  the  trough,  the  area  was 
still  en  the  back  side  of  the  trough  on  7  May.  A  new  upper 
level  trcugh  formed  and  approached  the  area  on  the  final  two 
days  cf  the  experiment. 

Surface  winds  were  relatively  light,  0  to  10  kts,  at  the 
beginning  cf  the  period,  and  increased  toward  the  end  of  the 
pericd  tc  16  kts,  gusting  to  22  kts. 

An  important  feature  of  these  interpretations  is  the 
nature  of  the  mixed-layer,  often  topped  by  an  inversion, 
with  reguard  to  stability,  and  therefore  mixing  intensities. 
It  is  assumed  that  the  mixing  becomes  greater  as  conditions 
become  mere  unstable. 
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At  the  start  cf  tie  experiment,  the  mixed  layer  exhibited 
stable  tc  slightly  unstable  conditions,  until  about  1800 
local  time  (PDT)  28  April,  when  conditions  became  mere 
neutral.  The  neutral  condition  remains  until  1  [lay  when 
conditions  ence  again  become  stable.  A  weak  frontal  passage 
before  0500  PDT  on  29  April  dees  not  appear  tc  affect  the 
nixed  layer  profile.  The  layer  remains  stable  until  a 
frontal  passage  on  2  May  when  conditions  become  neutral  and 
remains  so  until  5  May.  On  the  morning  of  5  May,  conditions 
are  slightly  stable,  but  returned  to  neutral  on  6  May, 
despite  a  frontal  passage  at  1300  PDT  on  5  May,  and 
remaining  neutral  to  the  end  of  the  experiment  on  9  May. 


C-   EiTA  SELECTION 

Data  frcm  the  MAGAT  experiment  were  selected  for  model 
verification  on  the  basis  of  the  relative  positions  of  the 
two  data  gathering  platforms  and  the  general  synoptic  condi- 
tions. The  locations  of  the  surface  ship  and  the  aircraft 
did  not  always  coincide  during  the  experiment.  Therefore 
soundings  and  aerosol  data  were  considered  suitable  if  the 
two  were  within  20  nautical  miles.  Hourly  wind  speed, 
direction,  and  SST  were  gathered  from  the  ACANIA  and  were 
used  alcng  with  the  aircraft  soundings.  Plots  cf  ' wind 
speed,  air  (solid  line)  and  sea-surface  (dashed  line)  temp- 
erature, and  relative  humidity  for  the  two  24  hour  periods 
are  shown  in  figures  3.9  and  3.10.  Another  consideration  of 
the  aircraft  data  was  the  relative  locations  of  the  aerosol 
ladders  and  the  soundings.  Although  both  gathered  during 
vertical  flight  profiles,  they  were  not  done  at  the  same 
time.  The  synoptic  conditions  during  the  experiment  include 
periods  cf  frontal  weather  which  was  also  avoided  in  the 
analysis.   Iwo  time  blocks 
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Figure   3.9        Plots  of  wind   spefd,    air  and   sea    surface 
temperature,    and   EH    beginning    1200   PDT   3   Say. 
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Same  as  figure  3.9,  except  beginning  1800,  6  May. 
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Figure   3.  11 


Flight   path   of  aircraft   on   the  morning   cf 
3   May,    1980. 
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Figure    3.12 


Same  as  figure  3.11  except  afternoon  of  3  May,  1980 
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Figure  3.13   Flight  path  of  aircraft  on  1  Bay  1980 
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Figure   3.14        Same  as  figure   3.13   except  6   Hay ,    1980 

39 


FLT#13 
7  MAY  80 


Figure    3. IS        Same   as  figure   3.13   except   morning    of   7    Say    1980. 
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Figure   3.16        Same   as  figure   3.13  except  afternoon   of   7   May    1980 

41 


+ 


+ 


MAGAT-80 


b  jaeaoeacd 


-f 


123c 


40' 


20' 


122' 


Figure    3.17        Cruise  track  of  R/V   ACANIA   3   May    1980. 
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Figure  3.18    Sai€  as  figure  3.16  except  H   Hay  1980. 
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Figure  3.19   Same  as  figure  3.16  except  6  and  7  Hay  1980, 
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were  chosen;  3-4  May  and  6-7  May.   Figures  3.11  through  3.16 

indicate  the   routing  of  the  aircraft   flights  for   the  two 

time  tlccks,  and  figures  3.17  through  3.19  depict  the  ccurse 

cf  the  R/V  ACANIA  during  the  same  times. 

At mcs pheric  data  from  the  aircraft  soundings  were  used 
instead  cf  the  radicscndes  frcm  the  R/V  ACANIA.  This  was 
done  fcr  twc  reasons:  first  because  the  aircraft  scundings 
contained  mere  levels,  and  second  they  were  in  closer  prcx- 
imity  tc  the  aircraft  aeroscl  ladders.  Atmospheric  and 
aeroscl  cata  for  the  first  run  was  taken  from  flight  6,  on 
the  ncrning  of  3  May.  The  locations  and  times  are  shewn  in 
figure  3.9.  The  aircraft  was  approximately  45  miles  frcm  the 
coast  when  a  sounding  was  taken  at  1143  PDT,  and  an  aerosol 
ladder  (111)  was  taken  at  1200  PDT.  Data  for  the  two  verifi- 
cations of  atmospheric  parameters  and  aerosol  values  were 
chosen  tased  on  their  proximity  to  the  initial  data.  A 
sounding  and  aeroscl  ladder  on  the  afternoon  of  3  May 
(sounding  time  of  17  14  PDT,  and  aerosol  ladder  113)  approxi- 
mately eight  hours  frcm  the  initial  time  was  chosen  for  the 
first  verif icaticn,  and  is  depicted  in  the  route  cf  flight 
7r  figure  3.10.  Data  for  the  second  verification  was  chosen 
from  flight  8,  figure  3.11;  sounding  time  1151  PDT,  and 
aeroscl  ladder  116.  Cata  for  the  two  verifications  are  bcth 
within  2C  miles  of  the  initial  data. 

Data  for  the  secend  model  run  is  shown  in  figures  3-12 
through  3-14.  The  model  was  initialized  with  data  frcm  the 
1732  PDT  scunding  and  ladder  122,  approximately  25  miles  eff 
the  Mcnterey  coast.  The  first  verification  data  was  taken 
from  flight  12  the  next  morning  approximately  17  hours 
later;  sounding  time  1743  PDT;and  ladder  123.  The  secend 
verification  data  was  taken  from  flight  13,  approximately  24 
hours  frcm  the  initial  data;  scunding  time  1836  PDT  and 
ladder  124.  Data  collected  for  the  two  verificatiens  were 
within  5  miles  of  the  location  of  the  initial  data. 
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Additional  wind  speed,  direction,  and  SST  collected  from 
the  B/V  ACANIA  were  used  for  initialization  and  verification 
cf  the  t*0  model  runs.  As  mentioned  in  the  model  descrip- 
tion, up  to  10  forecast  values  for  each  can  be  incut  into 
the  model.  However,  due  to  the  continuous  movement  of  the 
E/V  ACANIA  during  the  MAGAT  experiment,  where  SST  values 
were  changing  because  of  strong  coastal  gradients,  a  single 
value  of  SSI  was  used.  This  value  was  chosen  when  the  R/V 
ACANIA  was  closest  tc  the  areas  where  the  aircraft  data  was 
taken  for  each  cf  the  two  initialization  sites.  The  wind 
speed  and  direction  cata  from  the  R/V  ACANIA  was  considered 
representative  of  tiie  data  collection  sites,  and  therefore 
inpu*  into  the  model  runs  at  ac proximately  3  hour  intervals. 
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17.  MOD EI  INITIALIZATION  AND  RESULTS 

A.   DIGITIZED  ATMOSPEEBIC  SGUHEINGS 

Th€  atmospheric  inputs  for  the  model  (chapter  2)  ar<=: 

1.  The  nixed  layer  equivalent  potential  temperature  (C) 
and  specific  humidity  (gm/kg). 

2.  The  "jump"  discontinuities  in  each  of  the  above  values. 
The  junp  meaning  the  difference  between  the  mixed  layer 
value  and  the  value  at  the  top  of  the  inversion. 

3.  The  lapse  rate  for  each  parameter  value  above  the 
inversion. 

4.  The  depth  of  the  well  mixed-layer. 

To  simplify  the  process  of  calculating  these  values  from  the 
soundings,  a  digitizing  scheme  was  designed  by  the 
Environmental  Physics  Group  at  the  Naval  Postgraduate 
Schocl.  This  scheme  transforms  the  sounding  into  the  struc- 
ture that  describes  the  necessary  model  inputs.  Examples  are 
shown  in  figures  4.1  and  4.2.  Note  that  the  "jumps"  are 
depicted  as  occuring  in  an  infinitely  thin  layer. 

E.   AEBI1ICNAL  MODEL  IHPUT  AHE   ADJUSTMENTS  TO  THE  DIGITIZED 
SCUNEINGS 

The  reascn  for  making  adjustments  to  the  atmospheric  part 
cf  the  mcdel  is  to  provide  a  better  basis  for  evaluating  the 
behavior  of  the  aeroscl  prediction.  This  is  optimized  when 
the  mcdel  is  producing  the  best  possible  forecast  cf  the 
atmospheric  parameters.  The  parameters  are  forecast  very 
well  in  toth  model  runs,  with  these  adjustments  made.  The 
model  is  accurately  describing  the  height  of  the  inversion 
through  adjustments  in  the  subsidence  rate,  ana  the  lifting 
condensation  level   (LCL)   is  accurately  generating   clouds 
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through  adjustments  tc  the  temperature  and  specific  humidity 
which  are  within  the  limits  of  instrumentation  error.  The 
cloud  verifications  were  made  with  satellite  data  and  obser- 
vations frci  the  P./V  ACANIA  leg. 

An  initial  input  for  the  model  is  the  subsidence  rate. 
Using  a  first  guess  cf  -.0C5  ir/sec,  the  model  generated  a 
plot  cf  the  inversion  height  for  the  24-hour  period.  Eased 
on  the  verification  data,  the  subsidence  rate  was  further 
adjusted  tc  bring  both  the  predicted  and  observed  inversion 
height  values  together.  Then  a  small  adjustment  was  made  to 
the  mixed  layer  values  to  match  the  stability  and  cloud/ 
cloud  free  patterns  with  what  was  actually  observed  during 
the  period.  In  the  first  model  run  (3  May),  the  pctential 
temperature  was  increased  0.5  degrees  and  the  specific 
humidity  was  decreased  1  gm/kg.  The  resultant  sutsidence 
rate  necessary  tc  match  the  predicted  and  cbserved  inversion 
height  was  -.0042  m/sec.  These  adjustments  yield  an  atmos- 
pheric profile  similar  to  what  was  actually  observed.  The 
air  mass  became  slowly  saturated  enough  to  generate  clouds 
in  the  upper  part  of  the  mixed  layer  at  about  0300  PDT  on  4 
May.  Flcts  of  the  inversion  and  LCL  behavior,  and  the  temp- 
erature and  humidity  profiles  are  shown  in  figure  4.3.  In 
the  plot  for  LCL  and  inversion  height,  the  letter  "H" 
depicts  the  verification  times  and  values  for  the  inversion 
height,  and  the  letter  "L"  similarly  indicates  values  for 
the  LCL  verifications.  The  letters  "T"  and  "Q"  in  the  temp- 
erature and  specific  humidity  graphs  indicate  verification 
values  fcr  each  respective  value. 

The  atmospheric  plots  fcr  the  second  model  run  (6  May) 
are  seen  in  figure  4.4.  For  this  run  the  specific  humidity 
was  reduced  from  6.48  to  6.3  gm/kg,  the  potential  tempera- 
ture was  unchanged,  and  a  subsidence  rate  of  -0.01  m/sec  was 
input  for  the  run. 
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and  huaidity  profiles,    first   run. 
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Figure   4.4        Same    as  Figure    4.3,    except   for   second   run. 

52 


C.   AERCSCL  INITIALI24TIOS  EBOCEDURES 

The  aerosol  data  were  gathered  in  ladders,  usually  up  to 
five  fcilciceters ,  and  were  of  actual  concentrations  of  parti- 
cles for  a  given  increment  cf  radius,  n  (r)  .  This  number 
contained  bcth  marine  and  continental  components,  collected 
at  the  aibient  relative  humidity. 

The  initialization  scheme  was  performed  in  three  main 
steps.  The  first  step  was  to  identify  the  marine  and  conti- 
nental components.  By  using  the  digitized  soundings  for  each 
model  run,  a  deter  nination  was  made  of  two  height  values: 
cne  representing  the  well  mixed-layer  and  the  other  repre- 
senting the  free  troposphere,  above  the  inversion.  Eoth  the 
nixed  layer  and  free  troposphere  aerosol  counts  for  each  of 
the  five  radii  were  calculated.  The  second  step  was  to  apply 
the  humidity  growth  factor  to  the  values  and  calculate 
dv/dr .  The  relationship  between  n  (r) ,  which  was  measured, 
and  dv/dr  is  given  by  eguation  2-2.  The  final  step  was  to 
subtract,  for  each  rccius,  the  free  troposphere  dv/dr  value 
from  the  mixed  layer  dv/dr  value.  With  the  assumption  of  a 
constant  profile  of  continental  aerosols  through  both  the 
mixed  layer  and  into  the  free  troposphere,  this  subtraction 
leaves  in  only  the  mixed  layer  marine  components,  normalized 
to  80  percent  relative  humidity.  The  values  of  dv/dr  used 
to  initialize  the  model,  both  eguilibrium  and  MAGAT,  are 
shown  in  tables  4.  1  and  4.2. 

This  three  step  process  was  dene  for  both  model  runs,  and 
for  two  verification  times  within  each  run.  As  described  in 
chapter  III,  the  verification  ladders  were  chosen  to  be  as 
close  as  possible  to  the  position  of  the  initial  aerosol 
data.  In  both  cases,  the  verification  data  were  within  15 
niles  cf  the  initial  data.   The  plots  of  the  aerosol  output, 
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both  fcr  the  MAGAT  and  equilibrium  initializations,  is 
depicted  as  log  10  (dv/dr)  against  the  24  hour  time  period, 
in  half  hour  time  steps.  In  both  (initialization)  cases,  the 
"X"  values  represent  the  MAGAT  verification  values.  In  both 
model  runs,  aeroscl  output  for  the  MAGAT  initialization 
scheire  is  represented  by  a  continuous  line,  and  th€  dashed 
line    depicts  the  run    for   the   equilibrium    initialized   cases. 

TABLE    I 
Initial   equilibriui   aerosol   values 
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Wher  considering  the  equilibrium  initial  values  from 
table  one,  it  is  easy  tc  see  the  relationship  between  wind 
speed  and  production.  In  almost  all  cases,  the  aerosol 
volume  counts  increase  as  the  wind  speed  increases.  The 
distribution  of  radii  for  a  given  wind  speed,  however,  is 
not  as  well  defined.  For  lcwer  wind  speeds,  the  larger 
radii  are  few  ir  nunter,  and  as  the  wind  speed  increases 
there  is  a  dramatic  increase  in  the  number  of  larger  parti- 
cles. This  is  reasonarle  due  tc  white  cap  production.  Also, 
there  appears  to  be  a  slight  maximum  of  particles  at  2.0 
microns,  reguardless  of  wind  speed.  The  15  micron  initiali- 
zation values  in  the  tables  and  the  model  output  for  the 
same  si2e  are  not  shown.  This  is  because  the  observed 
initial  values  collected  from  the  aircraft  were  zero. 

In  comparing  the  MAGAT  initial  values,  one  finds  that  the 
smaller  radii  are  approximately  twice  the  equilibrium  count 
for  even  the  largest  wind  speed.  It  is  interesting  tc  note 
that  for  such  a  large  variation  in  initial  wind  speed  for 
the  two  days,  there  is  very  little  variation  in  the  initial 
values,  for  all  the  the  particle  sizes.  Based  on  results  of 
Monahan  /it  al  1983/  ,  it  was  also  surprising  to  find  no 
voluire  ccunts  for  the  larger  particle  sizes  since  wind  speed 
was  above  10  a/sec.  For  the  strong  wind  case,  there  were  no 
aeroscls  abcve  10  microns  collected  from  the  aircraft. 

D.   IIBSI  HCDEL  BON  BISOLTS 

The  flcts  of  dv/dr  over  the  24  hour  period  for  the  first 
model  run  are  shewn  in  figures  4.5  and  4.6.  For  each  radius, 
the  sclid  line  indicates  the  profile  with  the  MAGAT  initial- 
ized values,  and  the  dashed  line  represents  the  profile  with 
the  equilibrium  initialization. 
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For  the  three  sualler  sizes,  the  model  sho^s  little 
change  ever  the  pericd,.  however  the  model  does  producs  a 
trend  in  aerosol  production.  Whatever  volumes  were  used  to 
initialize,  the  number  tends  to  stay  within  one  crd«r  of 
magnitude  ever  the  period.  The  verification  values  cctrpare 
very  well  to  the  model  run  with  the  MAGAT  initialization. 
The  model  and  verification  dv/dr  values  differ  by  less  ^han 
10  units.  for  the  larger  size,  10  microns,  note  the  guick 
production  of  aeroscl  when  initialized  with  with  a  zero 
value.  Within  fcur  tcurs,  the  model  is  producing  clcse  tc 
the  same  values  as  the  eguilibrium  case.  From  this  fcur-hcur 
point,  the  two  cases  behave  very  much  the  same.  There  is  a 
diurnal  variation  evident  in  the  larger  size  cases;  a  lag  is 
noticeable  in  the  production  through  the  night  time,  with  as 
increase  after  sunrise  and  through  the  morning.  The  verifi- 
cation values  for  these  larger  sizes  are  not  good  for  either 
verification  scheme. 
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Figure  4.5        Aerosol   plots    for  first   run,    with  a)    r=  0.8   ui, 

and  b)     r=   2.0   um. 
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Figure   4.6 


Sane   as   figure   4.5,    except   a)    r~    5.0   urn, 
and  6)    r=    10   um. 
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E.       SECCKD    SODEL    EON    EESULTS 

Aerosol  plots  for  roth  the  equilibrium  and  MAGAT  initial- 
ization schemes  are  shown  in  figure  4.7  and  4.8.  The  dv/dr 
values  for  both  schemes  are  higher,  due  mainly  to  the 
increased  wind  speed.  For  all  of  the  par-icle  sizss,  there 
is  mere  cf  a  diurnal  variation.  This  is  slightly  evident  in 
the  smaller  radii,  and  well  defined  in  the  larger  radii. 
This  variation  is  strengthened  by  a  wind  speed  minimum  just 
before  sunrise.  Note  also  that  the  drop  in  production  of 
all  the  radii  around  C000  PDT  coincides  with  both  the  slight 
drop  in  wind  speed  frcm  10  to  8  m/sec,  and  the  formation  of 
clouds,      evident      in    the      LCL   and      inversion    figure.  This 

decrease  in  production  due  to  cloud  formation  is  not  notice- 
able in  the  first  model  run.  The  two  initialization  schemes 
are  ir  irore  agreement  than  with  the  first  run,  again  due 
mainly  tc  the  higher  wind.  With  the  larger  wind  speeds  in 
this  run,  the  model  is  generating  values  closer  to  the  veri- 
fication values  (in  tcth  initialization  schemes)  than  in  the 
first    run,    fchere   the    *ind   speed   is   significantly    less. 
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Figure   4. 8 


Same   as   figure   U.7#   except  a)    r=  5.0   ui   and 
b)     r=    10    um. 


61 


J.   SDMMiRI  OF  RESUI1S 

The  mcdel  is  producing,  with  a  minimum  of  adjustment,  an 
accurate  description  of  the  atmospheric  boundary  layer 
parameters.  Once  a  subsidence  rate  is  selected  such  that  the 
24-hcur  predicted  irversion  height  approaches  the  observed 
value,  the  values  of  temperature  and  humidity  only  require  a 
minimum  cf  adjustment.  In  bcth  cases  the  verification  values 
cf  temperature  were  within  2°  C,  and  the  values  of  specific 
humidity  within  1  gir/kg.  More  significantly,  the  trends 
ever  the  24-hour  periods  for  both  temperature  and  humidity 
were  accurately  described  by  the  model,  as  defined  by  the 
verification  data. 

With  respect  to  the  aerosol  input,  the  model  is  alsc 
producing  the  correct  trend  in  behavior.  There  is  evidence 
in  the  cutput  fcr  the  trend  in  aerosol  production  to  be  a 
function  of  not  enly  wind  speed  and  relative  humidity,  but 
also  the  concentration  of  the  initial  aerosol  values.  This 
is  seen  in  the  first  run,  particle  sizes  0.8  and  2.0 
micrcrs.  The  atmospheric  parameters  are  the  same  in  bcth 
cases,  fcr  each  radius;  but  it  is  the  initialization  scheme 
that  generates  a  different  trend  in  each  case.  The  mcdel  run 
with  the  MAGAT  initial  values  produces  a  loss  in  production, 
and  the  verification  values  are  within  one  crder  of  magni- 
tude cf  the  model  values.  With  the  equilibrium  case,  the 
trend  is  an  increase  in  production  over  the  time  pericd, 
with  the  mcdel  verification  values  two  orders  of  magnitude 
away  from  the  observed  verification  values. 

The  production  trends  in  the  second  model  run,  with  both 
initialization  schemes,  are  much  the  same.  In  this  case 
there  is  less  of  a  difference  between  the  model  verification 
values  fcr  the  two  initializations.  With  the  exception  of 
the  smallest  radius,   0.8   microns,   the  verification  values 
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for  all  the  radii  of  the  second  run  are  within  on?  crder  of 
aagnitude  cf  the  mcdel  output.  The  fact  that  the  MAGAT 
initial  values  are  closer  tc  an  equilibrium  value  at  a 
higher  wind  speed  seems  tc  imply  a  less  sensitive  mcdel, 
when  initialized  with  a  higher  wind  speed.  Accordingly, 
eeroscl  production  is  high;-  and  the  model  accurately 
predicts  the  behavicr  of  the  aerosol  in  the  stronger  wind 
case. 
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V.  CCHC IDS IONS 

Th€  boundary  layer  model,  with  a  minimum  cf  input,  will 
produce  an  accurate  description  of  the  MABL  [Davidson  ef  al, 
1984/  and  the  aeroscl  behavicr  within  the  marine  boundary 
layer.  !he  model  is  generating  values  for  marine  particles 
within  this  layer,  and  whenever  there  is  a  large  difference 
tetween  the  predicted  and  observed  values  ,  it  was  the 
observed  values  that  were  toe  high.  This  observation 
suggests  that  the  observed  values  collected  by  the  aircraft 
(made  up  of  marine  and  continental  particles)  may  have 
included  mere  continental  than  originally  expected.  A 
linear  distribution  (vertically,  through  the  boundary  layer) 
was  assumed  in  the  aerosol  initialization  process. 

Mere  significantly,  the  model  is  correctly  identifying 
the  production  trend  cf  aerosols  over  the  time  period.  The 
model  is  run  with  identical  atmospheric  inputs,  and 
different  aerosol  initialization  schemes  produce  two 
different  production  trends.  The  wind  speed  and  relative 
humidity  (the  only  atmospheric  inputs  of  previous  models) 
are  the  same,  and  the  model  is  generating  different  profiles 
based  on  initial  aercscl  concentrations.  With  a  minimum  cf 
two  case  studies,  it  is  clear  that  the  process  of  entrap- 
ment cf  air  from  above  the  mixed  layer,  containing  nc  irarine 
particles,  is  affecting  the  behavior  of  the  model  output. 
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